(N 



(N 



Sensitivity to Dark Energy candidates by searching for four-wave 
mixing of high-intensity lasers in the vacuum 

Kensuke Homma 1,2 

1 Graduate School of Science, Hiroshima University, 
Kagamiyama, Higashi- Hiroshima 739-8526, Japan 
2 International Center for Zetta-Exawatt Science and Technology, 



o: 

CN . Ecole Poly technique, Route de Saclay, Palaiseau, F-91128, France 

> 

q ; (Dated: November 13, 2012) 



Abstract 

Theoretical challenges to understand Dark Matter and Dark Energy suggest the existence of 



43 " low-mass and weakly coupling fields in the universe. The quasi-parallel photon-photon collision 

system (QPS) can provide chances to probe the resonant production of these light dark fields and 
the induced decay by the coherent nature of laser fields simultaneously. By focusing high-intensity 



lasers with different colors in the vacuum, new colors emerge as the signature of the interaction. 
Because four photons in the initial and final states interplay via the dark field exchange, this process 



> 

is analogous to four-wave mixing in quantum optics, where the frequency sum and difference among 



the incident three waves generate the fourth wave with a new frequency via the nonlinear property 
of crystals. The interaction rate of the four-wave mixing process has the cubic dependence on 
the intensity of each wave. Therefore, if high-intensity laser fields are given, the sensitivity to the 
weakly coupling of dark fields to photons rapidly increases over the wide mass range below sub-eV. 



X: 

Based on the experimentally measurable photon energies and the linear polarization states, we 
formulate the relation between the accessible mass-coupling domains and the high-intensity laser 
parameters, where the effects of the finite spectrum width of pulse lasers are taken into account. 
The expected sensitivity suggests that we have a potential to explore interactions at the Super- 
Planckian coupling strength in the sub-eV mass range, if the cutting-edge laser technologies are 
properly combined. 
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INTRODUCTION 



Ordinary matter occupies only ~ 4% of the total energy density of the universe. The 
remaining energies are supposed to be occupied by Dark Matter (DM) ~ 23% and Dark 
Energy (DE) ~ 73% In addition to the astronomical observations, directly probing 
these dark components in terrestrial laboratory experiments has crucial roles to provide 
different insights into the true characters of the universe or the structure of the vacuum. 

Light (pseudo) scalar fields are now indispensable theoretical tools to try to interpret the 

n 

cosmological constant A based on the DE scenario |2|. In reduced Planckian units with 
c = h = M P (= (8ttG)" 1/2 ~ 10 27 eV) = 1, the observed A ~ 1(T 120 is extremely small 
compared to the theoretically natural scale A ~ 1. There is a variety of theoretical models 
in the market. In order for a DE model to be falsifiable by laboratory experiments, we 
clarify following minimum requirements on the model: 

• solves the fine-tuning problem] how to realize such an extremely small A, 

• solves the coincidence problem] why the energy density coincides with the matter 
density only at once at present so accidentally among the long history of the universe, 

• predicts the field-matter coupling strength, the mass scale of the exchanged field, and 
the measurable dynamical effect, e.g., the force-range. 



For instance, quintessence approaches |2| are designed to resolve the fine-tuning and 
coincidence problems by introducing decaying behavior of A based on a potential of a scalar 
field. However, the potential forms are rather phenomenologically introduced. In the similar 
course, the scalar-tensor theory with A (STTA) Q], on the other hand, is grounded upon 
the conformal transformation and the frame of observations, which gives several testable 
predictions. Therefore, STTA is one of the DE models satisfying the above requirements 
simultaneously. The most significant prediction of STTA is the decaying behavior of A oc 
t~ 2 as a function of time t. The present age of the universe is to ~ 1-37 x 10 10 year 
corresponding to t ~ 10 60 ' 2 in reduced Planckian units. Thus, the observed A is naturally 



understood by the overal 



the dominant behavior The decaying behavior depends on the conformal frame on which 
our observations are based. For example, redshift measurements relevant to DE is implicitly 
based on the common atomic clock between distant points, i.e., on the common elementary 



decaying behavior, though expecting short-term fluctuations from 



particle masses. In order to realize constancy of particle masses, a consistent conformal frame 
must be favored on which the gravitational constant G looks constant, the expansion rate of 
the universe is consistent with the observation, and A decays as t~ 2 by keeping particle masses 
constant [7j. The choice of a conformal frame unavoidably associates a massless Nambu- 
Goldstone (NG) boson, because it breaks scale invariance (global conformal symmetry), 



which is also known as dilaton. In contrast to the well-known Brans-Dicke model jfj], a 
kind of scalar-tensor theories, the requirement of constancy of particle masses results in 
coupling of the scalar field with matter, i.e., violation of Weak Equivalence Principle only 
via quantum anomaly coupling [4j. Due to this coupling to matter fields, STTA uniquely 
predicts an extremely low-mass scalar field as a pseudo NG boson via the explicit symmetry 
breaking by the quantum loop effect in the self-energy. This is similar to a massive pion as 
a theoretical descendant of an originally massless NG pseudoscalar boson associated with 
chiral symmetry. The scalar field couples with other matter fields basically as weakly as 
gravity. The mass scale based on the simple one-loop diagram in which the light quarks 
and leptons with a typical mass m q ~ MeV couple to the scalar field with the gravitational 
coupling with the strength ~ Mp 1 is given by 

m 2 M 2 

m 4> J^2 ( 10 eV ) > (!) 

where the effective cutoff coming from the super-symmetry-breaking mass-scale M ss b ~ TeV 
is assumed, though allowing a latitude of several orders of magnitude, if M ss b is higher 



than the conventional TeV scale |4|. We note that the uncertainty on the mass range in 
the DE scenarios is quite large. Quintessence-based scenarios typically argue that the mass 
rom the second derivative of the assumed almost flat potential resulting in 



is determined 

m ~ 10~ 33 eV 46]. On the other hand, A ~ (meV) 4 in natural units intuitively leads models 
based on the particle picture jg- 10] assuming the mass scale in the meV range via rather 



complicated assumptions. For example, the axion inspired models [8] share the similarity to 
STTA by introducing the concept of pseudo NG boson driven by the two dominant scales, 
Mp and a scale of the assumed symmetry breaking at a lower energy than Mp |47l| . 

The finite mass of the scalar field in STTA causes non-Newtonian force [11] via the 
Yukawa potential, a.k.a. fifth force. The inverse of Eq.flTJ gives a finite range corresponding 
to ~ 100 m of the force mediated by the exchange of a quantum between local objects. 
This is an entirely different aspect from its way of a cosmological fluid in accelerating the 



universe. The force-range has not been explicit in the other theoretical models. This unique 
aspect triggered the past experimental efforts to measure deviations from the Newtonian 
potential between massive test bodies 12| in different contexts from DE at that time. These 
measurements, however, accompany large systematic uncertainties due to the uncontrollable 
macroscopic probes. As an alternative approach, we have proposed to utilize the nature of 
high-intensity laser fields toward laboratory search for the scalar field predicted by STTA [k^ 
as an ultimate goal of laboratory experiments. 

Furthermore, low-mass and weakly coupling fields are also predicted in the contexts of 
particle physics with the solid foundation. For example, axion, the pseudoscalar field is pro- 
posed as a NG boson associated with the global Peccei and Quinn symmetry breaking [14J] 
to naturally maintain the CP conserving nature of the QCD Lagrangian. Axion and invis- 
ible axion-like fields have been intensively investigated by astrophysical objects as well as 
laboratory experiments jl5|. Some of them may become cold dark matter candidates, if the 
mass and the coupling to matter are within the proper range. Such fields may also leave 
observational isocurvature fluctuations, if the symmetry breaking occurs during the infla- 
tion phase of the early universe 16j. If we can anticipate that the experimental sensitivity 
reaches the gravitational coupling strength, the detection of such cold matter candidates 
with much stronger couplings to matter naturally comes into view as the preliminary step 
toward the ultimate goal. 

We, therefore, generalized the principle of the measurement to search for both scalar and 
pseudoscalar fields in a model independent way as much as possible [17T | . As amplitudes 
of laser fields are increased, we can improve the sensitivity to weakly coupling low-mass 
fields predicted by any types of theories, as long as the coupling to photons is expected. 
The proposed method can be regarded as a kind of particle colliders attempting to produce 
extremely light resonance states. The mass range of interest is, however, for instance, much 
lower than that of Higgs-like boson produced at Large Hadron Collider (LHC) by more than 
ten orders of magnitude. In the proposed method, following two key ingredients to enhance 
the sensitivity are included. 

The first ingredient is the introduction of the quasi-parallel photon-photon collision sys- 
tem (QPS) as illustrated in Fig{TJ This is considered to realize the center-of-mass system 
(CMS) energy as low as possible between colliding two laser photons for the production of a 
low-mass field as a resonance state without lowering the incident photon energy below opti- 
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cal frequency. The head-on collision in CMS corresponds to the case i? = 7r/2 in Figj3j The 
quasi-parallel system can be obtained by introducing a Lorentz boost of the head-on collision 
into the perpendicular direction with respect to the incident direction in CMS. The CMS 
energy E cms is then expressed as E cms = 2a;sin'# with the incident energy of photon u. If a 
small t? is realized in the laboratory frame, it provides an extremely low CMS energy. This is 
the essence of the introduction of QPS. Moreover, fortunately thanks to the strong Lorentz 
boost, QPS provides frequency-shifted photons as the decay product of the resonance state, 
which becomes a distinct observational signal as the indication of the interaction. However, 
as we discuss in the next section briefly, the resonance point cannot be directly captured due 
to the extremely narrow resonance by the weakly coupling compared with the momentum 
uncertainty of incident photons in QPS. This situation requires an averaging process of the 
cross section over the possible uncertainty of the CMS energy in QPS. By this averaging, 
the non-negligible effect of the narrow resonance is enhanced by the square of the inverse 
coupling compared with the case where no resonance state is contained. If the coupling to 
two photons is proportional to 1/Mp, the huge enhancement by Mp is expected. 



The second ingredient is the enhancement by the coherent nature of laser fields 19[ or 
the degenerate nature of Bosonic particles as illustrated in FigfSJ This Bosonic nature of the 
laser beam is fundamentally important, because we can induce the decay of the produced 
resonance state into a specific momentum space as the principle of the laser amplification 
itself utilizes that nature. We propose to use different frequencies between the production 
and inducing laser beams, respectively. As shown in the figure, the exchange of the low- 
mass field is interpreted as the four-wave mixing process where three waves (the two waves 
are degenerate and the one wave has a different frequency from the degenerate waves) are 
combined and the forth wave emerges with a new frequency not included in the origi nally 
mixed laser waves. This four- wave mixing process is well-known in quantum optics [20|. 
The process is already applied to generate a different color wave from those of incident laser 
beams via the nonlinear atomic processes of crystals. In other words, the proposed method 
is as if the atomic nonlinear process is replaced by the nonlinear process of the vacuum via 
the low-mass field exchange. In the context of the QED interaction, a similar approach is 
discussed [2l| and the experimental setups are proposed [22| . The upper limit of the photon- 



photon cross section is provided by this method 



231 ] . Since each of the two photons at the 



first vertex annihilates into the coherent state with lu, while another photon at the second 
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FIG. 1: Quasi parallel colliding system (QPS) between two incident photons out of a focused 
laser beam with the focal length /, the beam diameter d, and the upper range of incident angles 
A$ determined by geometric optics. The signature (2 — u)uj is produced via the four-wave mixing 
process, lu + lu — > (2 — u)uj + ulo with < u < 1 by mixing two waves with different frequencies 
Iuj and uuj in advance at the incidence. 



vertex is created from the coherent state with uoo with < u < 1, the interaction rate to 
observe (2 — u)oj frequency is eventually enhanced by a factor of (y/ N luJ y/ Ni w ^ N UUJ ) 2 [131 ]. 
where N indicates the average number of laser photons with individual frequency specified 
by the subscripts. This cubic dependence of the interaction rate motivates us to make the 
laser energy per pulse as large as possible. 

nn 

In this paper, we extend the formula discussed in the recent works [13J, H7| and then 
provide the prescription to relate the accessible mass-coupling domains by taking an essence 
of realistic experimental constraint such as a state of the multi-frequency mode with a finite 
frequency bandwidth and the effect of the specification of linear polarization states, when we 
attempt to apply this method to experiments based on pulse lasers. The expected sensitivity 
to the low-mass and weakly coupling fields is provided for anticipated high-inte nsity laser 
fields available at laboratories over the world at present and in the near future 24j. The 
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FIG. 2: Four- wave mixing process in the vacuum. The state \N)) = e N '^_Y^L 



N n/2 



n) with 



19| with the average 



\n) being the degenerate state of n photons refers to the coherent states 
number of photons N . The coherent states are distinguished with photon frequencies specified by 
the individual subscripts. The probability to create (2 — u)uj with < u < 1 is expected to have 
the cubic dependence on N, because expectation values of creation and annihilation operators 



associated with the individual photon legs become \/N in the scattering amplitude [13]. The 
effective dark field - two photon coupling gM~ l for each vertex corresponds to the coupling in the 
effective interaction Lagrangian defined in Eq.([2|). 



result suggests that the state-of-the-art technology may provide access to interactions with 
gravitational coupling strength and even beyond it (Super-Planckian coupling) for relatively 
higher mass ranges in the sub-eV mass domain. We emphasize that the proposed approach 
is a kind of Bosonic collider. The commonality and the distinctions from the Fermionic 
collider, for example LHC, is discussed as a concluding remark. 
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FORMULAE TO RELATE SENSITIVITY AND LASER PARAMETERS 



Let us briefly review the necessary parts for the extension in order to consider the effects 
of the finite spectrum widths and the specification of linear polarization states of laser fields. 

The effective interaction Lagrangian between two photons and an unknown low-mass 
scalar or pseudoscalar fields or a are generalized as follows, respectively 

-L^gM-^F^y or - L a = gM^-F^a, (2) 

where M has the dimension of mass while g being a dimensionless constant. Depending on 
the allowed polarization combinations of two photons coupling to the dark fields, we can 
argue whether they are scalar-type or pseudoscalar-type in general as we see in Appendix 
in detail. 



We summarize the notations and kinematics based on the equation (2.1)-(2.3) of 13]. We 
label momenta to four photons as illustrated in Figj3j where the incident angle i9 is assumed 
to be symmetric around the z-axis, because we assume the symmetric focusing as illustrated 
in FigJU For later convenience, we introduce an arbitrary number u with < u < 1 to 
re-define momenta of the final state photons as = uuj and 0J3 = (2 — u)u. By this 
definition, we require < W4 < U3 < 2u>. We consider the case where we measure U3 with 
the specified polarization state as the signature of the interaction. With these definitions, 
energy-momentum conservation in [l3j is re-expressed as 

(2 — u)u) + uuj = 2uj (3) 

(2 — it)u; cos #3 + mo cos 6± — 2u cost? (4) 

(2 — u)u sin 9s — uuj sin #4. (5) 
From these, we also derive the following relation 



sin 63 sin 2 $ 



sin 64 1 — 2 cos $ cos 64 + cos 2 d ' 

and 



(6) 



^sin 2 ^ 

w 3 = (2-«a;=- -. 7 

1 — COS V COS 63 
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From (2.5) of 13], the differential cross section per solid angle of p% is expressed as 

,2 /,.,_\2 

6 3 \87rwy 

where |A^5| 2 is the square of the invariant scattering amplitude including the resonance state 
in the s-channel with a sequence of four-photon polarization states S = Pi^Ps,^ specified 
by the polarization vectors with i = 1,2,3,4 for the photon labels, whereas — 1,2 
are the kind of the linear polarization as indicated in FigfJJ As we discuss how to evaluate 
the amplitudes in Appendix, S = 1111, 2222, 1122, 2211 and S = 1212, 1221, 2112, 2121 give 
the non-vanishing invariant amplitudes for scalar and pseudoscalar exchanges, respectively. 

The enhancement by the inducing laser field labeled as p^ is limited to the intrinsic 
spectrum width of the inducing laser energy due to energy-momentum conservation, which 
is defined as 

5u4 = ull — uji = (6u)u), (9) 

where 5u = u—u with 0<m<1,0<m<1 and u > u. The overline and underlines attached 
to the simple variables indicate the allowed maximum and minimum values, respectively. 
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This notation is used repeatedly, unless confusion occurs. The spectrum width of u 3 is 
simultaneously constrained to 

5u 3 = uj 3 — uj 3 = (2 — u)u — (2 — u)oj = (5u)uj (10) 

due to energy-momentum conservation as well. This indicates that the energy range to 
detect Co> 3 must be consistent with <5w 4 in order to get the full enhancement factor by the 
inducing laser. 

By assuming that 5u 3 = 6004 is realized in an experimental setup, let us now calculate 
the cross section integrated over 5u> 3 as a function of the spectrum width parameter (u,u). 
Eventually 5u) 3 gives the range of the integral via the range of scattering angle 9 3 . The 
integrated cross section over the range of 9 3 from 9 3 to 9 3 is expressed as 



Fs\M s \ 2 [°* (u 3 \ 2 . a JO _ Fs\M 



2 



*= ,0 / 40 ^ S m9 3 d9 3 = X, (11) 



(12) 



(8W) 2 sm 4 tf Je, \2ujJ 00 (8tto;) 2 

by defining 

1 fa /cu 3 \ 2 r 1 183 

X=—^ / f-i] sm9 3 d9 3 



sin 1? J 3 \2uJ 



4 cos #(1 — cos 9 3 cos •&) 



s 



where |.Ms| 2 denotes the averaged |A^s| 2 over possible uncertainty on the incident angle •& 
as we briefly discuss below, and JFs corresponds to the integral over the azimuthal degree of 
freedom depending on the specification of photon polarizations in the initial and final states. 
If the scattering amplitude has the axial symmetry around the z-axis in Fig|3j J-$ simply 
corresponds to 2n. As we summarize in Appendix, however, the axial asymmetries actually 
appear depending on S specified by experimental conditions, which result in deviations from 
2tt. We then convert the variable cos #3 to u based on Eq.fJZJ) from which we express 

cos fl 3 = _J_ (1 _ ^ s inH) ~ 1 + \& (l - —) = 1 + U 2 (l - —) , (13) 



cos $ V cu 3 J 2 \ tu 3 J 2 V 2 — u, 

where terms of the order higher than d 2 are dropped when applied to the low-mass resonance. 
From this, Eq. (112j) is also approximated as 

X ~ — (14) 

We now consider the scattering amplitudes only for the case when low-mass fields are 
exchanged via the resonance states in the s-channel. The resonance decay rate of the low- 
mass field with the mass m into two photons is expressed as 

T = (IQn)- 1 (gAr 1 ) 2 m 3 . (15) 
10 
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As we calculate in detail in Appendix, for example, in the case of scalar field exchange, 
the invariant amplitude in the coplanar condition where the plane determined by p\ and P2 
coincides with that determined by p 3 and p 4 is expressed as 

/ ,^ is9 ^ 4 (cos 2$ — l) 2 
Ms = -(gM ) - \ -j- 2 , 16 
2u z (cos M — 1) + m z 

where the denominator, denoted by T> in the following, is the low-mass field propagator. We 
then introduce the imaginary part due to the resonance state by the following replacement 

m 2 — > (m — iV) 2 « m 2 — 2imT. (17) 

Substituting this into the denominator in Eq. ( |T6|) and expanding around m, we obtain 

V » -2 (1 - cos 2#) (x + io) , with x = w 2 - w 2 , (18) 

where 

2 m 2 /2 mT 
Ur= l-cos2tf' a= l-cos2^- ( } 

From Eq. ( TT5|) and (fT9j) . a is also expressed as 



W£ /gfflV 
8tt UJ 

which explicitly shows the proportionality to M~ 2 . We then express the squared amplitude 



a=^r—\ , (20) 



as 

,2 



l-^sl 2 « ( 47r ) 2 7TT^- ( 21 ) 



a 

X z + a" 

Theoretically if we take the limit of w — > u r , |A^5| 2 — > (47r) 2 is realized from (12TI) . This is 
independent of the smallness of a cx M -2 . Meanwhile, the off-resonance case \ ^> a, |^Vf ,5 1 2 
equivalent to Eq.( jT6i) is largely suppressed due to the factor a 2 oc M~ 4 for the case of a 
weakly coupling M -1 . 

This is the most important feature arising from the resonance that overcomes the weakly 
coupling stemming from the large relevant mass scale such as M = Mp. However, we are 
then confronted with an extremely narrow width a for e.g. gm 1 eV , M ~ Mp = 
10 27 eV and u r ~ 1 eV. The rescue to overcome this difficulty is the averaging process over 
unavoidable uncertainties of incident angles in QPS. Even if a single photon with a fixed 
frequency u is focused by a lens element in QPS, the wave vector around the diffraction 
limit fluctuates by the wavy nature, in other words, the beam waist Ax at the diffraction 
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limit and the momentum accuracy Ap must satisfy the uncertainty principle AxAp > h/2. 
Therefore, we need to distinguish between the theoretically specified momenta and the 
physically specified ones in QPS. 



Based on Eq.(40),(41), and (42) of [17j], we express the average of the square of the 
invariant amplitude over a possible uncertainty on the incident angle i.e., uncertainty on 
the directions of the wave vectors as 

]Ms? = M)\Ms\ 2 d* ~ ^ &y«, (22) 

where the resonance angle $ r satisfies the resonance condition u 2 = u r = (m 2 / 2)7(1 — 



cos2i? r ), and we plugged the following simplest angular distribution function p [48j] into 



(23) 



1/Ad forO<$<A?9 
for Atf < i? < tt/2 

which is normalized to the physically possible range < d < tt/2. The incident angle 
uncertainty A$ can be as large as that determined from geometrical optics j^] 

A?? ~ ^, (24) 

with the common beam diameter d and focal length / for both the creation and inducing 
beams as illustrated in FigJTJ 

By substituting Eq.( IT4"]) . ( )22|) and (!20|) into Eq.( ITT|) . we express the partially integrated 
averaged cross section as 



__ J r s\M s \ 2 Fs (4tt)2 uP_ /gm\ 2 bu_ 

" ^ ~ ' 2u 2 \A$ IQtt\m) n '4$ 2 



-5 bu f gm\ z 1 / A \ 2 bu I gm\ 
512u 2 AA^ \ MJ ~^Y2\2^j tirAfi S \M J 



Ts Su ( qm\ 2 1 ( A \ bu ( ' gm\ 2 



In the averaging process, among a possible range of only d ~ d r effectively contributes 
to the cross section because of the narrow width a in Eq.( l2"Uj) for a large M. The second 
line in Eq. (125|) takes this aspect into account. For the last equation, w[eV] = 2irhc/X with 
h = c = 1 is substituted, where A is the wavelength of the creation laser field. 

We now express the yield of frequency-shifted photons y as a function of the spectrum 
width parameter (u, u) 

y = Ca (26) 
12 



with the effective integrated luminosity £ over propagation time of a single shot laser fields 
with pulse duration time r which is assumed to be common for both the creation and induc- 
ing beams. We discussed how the integrated effective luminosity should be defined in 
Here we briefly review the relevant part. The solution of the electromagnetic field propaga- 
tion in the vacuum with a Gaussian profile in the transverse plane is well-known |l8[]. The 
transverse spatial profile of a laser field is typically Gaussian to the first order approxima- 
tion in high-intensity laser systems. In this case, the electric field propagating along the 
z-direction in spatial coordinates (x, y, z) is expressed as 

E(x,y,z) oc 

wM-' lkz - Hiz)] -'i^ + ^m)}- (27) 

where k = 2ir / X, r = y/x 2 + y' 2 , wq is the minimum waist, which cannot be smaller than A 
due to the diffraction limit, and other definitions follows: 

w{zf = w 2 (l + , (28) 



R = z(l + Z -f), (29) 



H{z) = tan" 1 — , (30) 
\z R J 

_ ™ 2 

zr = (31) 

The transverse beam size of the focused Gaussian laser beam is minimized at the beam waist 
and then expands beyond the focal point where interactions among out-going photons are 
prohibited by the condition that two photons propagate into opposite directions in CMS. 
On the other hand, the exchange of a low-mass field may take place anywhere within the 
volume defined by the transverse area of the Gaussian laser times the focal length / before 
reaching the focal point (see FigJT]). 

Given the Gaussian laser parameters above, the effective integrated luminosity C over 
the propagation volume of a laser pulse can be defined as follows \y\ . At an instant, the 
interaction is limited within a region over cr where the average number of photons iV c and 
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Ni are available for creation and inducing processes, respectively. Luminosity at a point z 
integrated over pulse duration r is expressed as 

m = Mg>> = ^4 (32) 

where I(N C , N) denotes a dimensionless intensity depending dominantly on the average 
number of creation and inducing photons, N c and Ni, respectively within duration time 
r. The expression w 2 (z) in Eq. ( 1281) is substituted. During the propagation over the focal 
length /, the effective number of the interacting regions or the number of virtual bunches 
b is expressed as b = //(cr). Therefore, the effective integrated luminosity C over pulse 
propagation time averaged over the focal length / is finally expressed as 

In the case of charged particles or Fermionic beams, the dimensionless intensity I corre- 
sponds to the square of the number of particles per bunch, which is the combinatorics to 
take two Fermions from individual colliding beam bunches. On the other hand, in the case of 
four-wave mixing, all photons are in the quantum coherent states with the inducing nature 
resulting in the cubic dependence as we discussed. By taking this aspect into account, we 
define the dimensionless intensity included in C as follows 

UNcNJ^CrriNtANi, (34) 

where C m b is a factor to consider combinatorics for the choice of two photons in the creation 
beam and one photon in the inducing beam by extending the argument for the single- 
frequency mode [13!] to the multi-frequency mode as discussed below, and A is an acceptance 
factor for the inducing photons to satisfy energy-momentum conservation in the final state. 

The acceptance factor A is introduced because the process occurs only in a small portion 
504 of the entire angular spectrum with the whole strength N distributed over the total 
range of A6* 4 , hence; 

A= Af, = M> (35) 

which is much smaller than unity. Here A# 4 is further assumed to be common with that of 
the creation beam A$ by sharing the same optics as that of the creation beam. The 56$ is 
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constrained by the spectrum width of 0J4 = uu, therefore, is described as a function of the 
spectrum width parameter (u, u) as follows 

This relation is obtained from energy-momentum conservation as follows. Equation (jSJ) gives 
the ratio 7Z defined by 

^ sin #3 oja u ,„„. 

71 = / = — = . (37) 

sin 64 U3 2 — u 

By equating Eq. (j3"7j) and (jSj), we obtain 

ft(l + cos 2 $) -sin 2 tf , . 

cos# 4 = — — '— . 38 

2ftcostf v ; 

Neglecting higher order terms more than -# 2 , we approximate 84 as 



<WW^- (39) 

For a low- mass case m ~ 2'doj with •& -C 1, Q4 is also small via 64 ~ 7lr l ^ 2r d from Eq. ()39l) . 
Eventually the emission angle of the signal, #3 also becomes small via 9 3 ~ ft 1 / 2 ?? from 
Eq. ( 13 7p and fl39|) . This is the reason why the creation and inducing beams are all aligned 
into the same optical axis z as illustrated in FigJU by which a chance to enhance the inducing 
process is maximized. We note here that we have only to search for (2 — u)u as the signature 
of the four-wave mixing process without measuring the emission angle #3 directly. 

We then consider the effect of multi-frequency mode of a short pulse laser via the fol- 
lowing argument on the combinatorics factor C m b in the mixing process. This is physically 
unavoidable, because a pulse laser with a short time duration must contain the correspond- 
ing energy uncertainty in principle. We assume that photons of a creation laser pulse in 
a state of multi-frequency mode can be uniformly divided into n c frequency bins, namely, 
forming an uniform frequency density within the frequency bandwidth of the creation laser 
pulse, each of which forms a coherent state of the individual frequency. We also assume 
rii for the inducing laser pulse as well. All frequencies as a result of the possible mixing 
between creation and inducing frequencies must be detectable by an experimental setup 
in order for the following argument to be valid. We also assume all the frequency modes 
share the common focusing length and beam diameter independent of the wavelengths as 
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illustrated in FigJTJ In this model case, the dimensionless intensity in Eq.( l34"J) included in 
the luminosity definition is expected to be 

1 9 / [n~ c [n~ c \ 2 ( /iv7\ 2 1 9 

2 V V n c\ n c J \\j rii J 2 

where square roots are enhancement factors due to the coherent states of individual frequency 
bins, \n 2 c corresponds to combinatorics to choose initial two photons within n c frequency bins 
of a creation laser pulse by allowing to choose two photons even within the same frequency 
bin, while rii is the degree of freedom to chose one frequency bin out of rii bins. We thus 



find C m b = \ which is the same dimensionless intensity as discussed in 17| after all. 

The effect of the multi-frequency mode is not seen directly in terms of the combinatorics, 
however, we put a note that choosing two photons uj\ and U2 out of n c bins implies that the 
incident two photon energies could be different, which deviates from the simplest assump- 
tion of the symmetric angle of incidence with equal photon energies as illustrated in FigJTJ 
However, we can always find a reference frame by a Lorentz boost which exactly satisfies 
the symmetric conditions by considering the inverse process of the massive particle decay 
into two photons starting from the rest frame of that particle, because the interaction is 
enhanced only when the resonance condition m ~ 2-du is fulfilled. Instead, p% in QPS must 
contain fluctuations by this boost effect. Therefore, as long as the experimental coverage on 
p 3 with respect to the additional fluctuations is broad enough, the integrated interaction rate 
over the detector acceptance is unchanged from the case with single mode lasers. In order 
to estimate the necessary coverage of u 3 approximately, we may re-define uj = (u>i + w 2 )/2 as 
the averaged value between arbitrarily selected two incident photon energies by taking the 
approximation where the transverse momentum of the produced massive particle is negligi- 
bly small compared to the longitudinal momentum of the produced particle in QPS, and all 
equations in Eq.([3])-([6]) are restored, because they are simply scaled by the newly defined oj. 
Therefore, the arguments so far are approximately valid even in the case of the multi-mode 
lasers. Suppose that the bandwidth of the creation beam and inducing beams are defined, 
respectively as 

(u)-Aw<w< (u) + Au (41) 
(w 4 ) - Aw 4 < u 4 < (0J4) + Aui, (42) 
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where ( ) denotes an average of each frequency distribution, while A indicates half of the 
full bandwidth of each frequency. According to the last part of Eq.(JH]), Su^ varies depending 
on an arbitrary chosen u. However, u is actually defined with respect to u. Therefore, 
any chosen u are absorbed into the definition of u = uj^/uj. Thus only the first part of 
Eq.(JHD becomes essentially relevant, which is determined by the intrinsic character of the 
prescribed inducing laser beam independent of the creation laser frequency Ioj. Accordingly 
an experiment should be designed so that all the modified range of 5u 3 via relation Eq.(j9]) 
and (flOl) is acceptable. Let us introduce new notations to describe the modified range of 
o>3 due to the shift of the averaged u in order to distinguish it from 5u 3 intrinsically caused 
by 5w 4 via Eq.flU]) and (IT0|) . The modified upper and lower edges of u 3 can be defined, 
respectively as 

u 3 = 2W - W4 = 2((u) + Au) - ((w 4 ) - Aw 4 ) = (co 3 ) + Au 3 (43) 

^3 = 2a; - ull = 2((lu) - Aw) - ((w 4 ) + Aw 4 ) = (u 3 ) - Aw 3 , (44) 

with (u 3 ) = 2(u) — (o; 4 ) and Au 3 = 2Au + Ao; 4 . As long as an experiment can accept u 3 in 
the range u' 3 < lu 3 < lu' 3 , the expected yield in the proceeding paragraphs is valid. Naturally, 
the range of Soj 3 via the intrinsic Suj^ is fully contained in this modified range. In the case of 
the multi-frequency mode, strictly speaking, u should be defined as u = (w 4 ) / (u>) resulting 
in u = ((w 4 ) + Acj 4 )/ (a;) and u = ((w 4 ) — Aw 4 )/ (u). For the entire arguments in this paper, 
u and 5u are implicitly assumed to be defined based on the averaged frequencies. 
We now express the yield y by substituting Eq. flM]) . fl35|) and ff25|) into Eq.f l26|) 

y = K (X, r)Kx(f, d)K 2 (u,u)J : 's (j^j ' C mb N c 2 N h (45) 

where 

K (X,t)^- (46) 
cr 

indicates that shorter pulse duration time cr — > A has the maximum gain on the yield, 

™ - soisW-" (i) ~ «b (ff^ 1 (i) (47) 



with wo = f X/ (-n(d/2)) for zr = tiwq/X |l8| is the parameter relevant to only optics, and 




K 2 (u,u)=USu= [J =-J^—\(u-u) (48) 
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is the laser spectrum width parameter determined from the spectrum width of the inducing 
laser field, 5u4. 

From Eq. fl45p we finally obtain the expression for the coupling parameter g/M to discuss 
the sensitivity as a function of m for a given experimental parameters via the following 
equation, 



For convenience to design experiments, we also consider the case where duration times are 
not equal between the creation beam r c and the inducing beam Tj. Since the yield is increased 
by the quadratic dependence on the creation beam intensity, it is natural to realize the case 
t = t c < Ti in experiments. If this is the case, the accessible coupling is re-expressed as 



EXPECTED SENSITIVITY 

First we summarize the key control parameters or experimentally adjustable knobs based 
on the arguments in the previous sections. The resonance condition is satisfied if m ~ 2du. 
However, instead of hitting the resonance point directly, our approach is to take the average 
of the squared scattering amplitude over the possible incident angle uncertainty A$ in QPS 
by the focused laser beams. Changing the focal length introduces different A$, i.e., the 
different range of the angular integral for the averaging. If a resonance peak is contained in 
that range, the resonance effect appears as the integrated result. The basic strategy of this 
proposal is therefore to change the focal length, attempting to search for the appearance of 
four-wave mixing photons, which approximately gives a mass range via the relation m < 
2uAi}. Thus, this knob provides variations along the mass axis, m. On the other hand, the 
four-wave mixing yield is enhanced by the cubic product of the laser intensities. Therefore, 
changing laser intensities gives large variations along the coupling axis, g/M. If a significant 
signature is found, we can localize the domain in the (m, g/M) plane by adjusting these two 
knobs. 

In addition, there are knobs on the polarization states of laser fields. As we discuss 
in detail in Appendix, depending on the types of exchanged fields, different polarization 
correlations are expected between the two photons in the initial and final states. In the case 




(49) 




(50) 
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of the scalar field exchange which is the first of Eq.([2l), the possible linear polarization states 
in the four-wave mixing process are expressed as follows: 



co{l} + tu{l} -)• (2 - u)lo{1} + uco{l} 
+ ->■ (2 - u)u{2} + uu{2}, 



(51) 



where photon energies from the initial to final state are denoted with the linear polarization 
states {1} and {2} which are orthogonal each other. On the other hand, in the case of the 
pseudoscalar field exchange with the second of Eq.(J2]), the possible linear polarization states 
are expressed as: 



co{l} + co{2} ->• (2 - u)u{2} + ulo{1} 
+ u{2} ->■ (2 - u)u{l} + uu{2}. 



(52) 



By choosing physically allowed combinations of the linear polarizations, we can distinguish 
the types of exchanged fields, while we can estimate the background processes by requiring 
the false combinations on purpose. 

Let us now briefly review some of major high-intensity laser facilities in the world in- 
cluding on-going projects which can provide more than 100 J per pulse. There are typically 
two classes of laser systems to achieve high-intensity: a moderate pulse energy per tens of 
fs short duration and a large pulse energy per several ns duration. The established choices 
of the laser technology for the former and latter classes are Titan:sapphire-based lasers and 
Nd:glass-based lasers typically dedicated for laser fusion studies, respectively. We note that 
energy per shot is more important than pulse duration for the proposed method, because 
the interaction rate is cubic to the numbers of photons, while it is inversely proportional to 
pulse duration. On the other hand, the typical repetition rates for such high-energy pulse 
lasers are currently limited to at most every minute and every several hours for the for- 
mer and latter classes, respectively. The frontiers of the former class are VULCAN 10PW 
300J/30fs[32|, APOLLON 150J/15fs[33|, and what is prepared for the Extreme Light In- 
frastructure (ELI) project by combining the 20 APOLLON-type lasers 1 34]. The examples 
of the latter class are GEKIKO-XII 0.1- lkJ combining 12 beam lines [sj, VULCAN 2. 6k J 
combining 8 beam lines 36[, FELIX 10k J combining 4 beam lines 35], OMEGA 30k J com- 



bining 60 beam lines 37|, PETAL with quad LMJ 1.5kJ-80kJ 



240 beam lines 



39], and NIF 1.8MJ combining 192 beam lines 40|. The quoted numbers 



LMJ 1.8MJ combining 
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above should be regarded as rough references which can change depending on the opera- 
tional conditions and the progress of the on-going projects. The bottleneck of the currently 
available laser systems with respect to the proposed method, the higher the pulse energy, 
the lower the repetition rate, is going to be improved. The proposed technique of Coher- 
ent Amplification Network (CAN) [41] adopts the coherent addition of highly efficient fiber 
lasers which is in principle operational at a higher repetition rate resolving the heat problem 
typically seen in the Nd:glass-based lasers. An international project, International Coherent 
Amplification Network (ICAN), has been launched 421 ] . This is also encouraging for the 
community of high energy physics aiming at much higher energy than that of the conven- 
tional acceleration technique, because the high average power is eventually necessary for high 
luminosit y p hysics even based on the new acceleration scheme; Laser Wakefield Acceleration 
(LWFA) [43]. High-intensity lasers can serve for the extension of the traditional course of 
high energy physics as well as the novel type of physics as discussed in this paper. Interna- 
tional Center for Zetta- and Exawatt Science and Technology (IZEST) has been launched 
toward the integration of high energy physics and high-field science [45|. These combined 
efforts are reviewed in our recent article 44]. 

Figure H] indicates explorable domains in the (m,g/M) plane by searching for the four- 
wave mixing process by counting the number of photons in the frequency-band u' 3 < uo 3 < u' 3 , 
where J-$ = 2tt is used because any J-g are on the same order as those in the axial symmetric 
case as we discuss in detail in Appendix and the intention of this plot is not in the separation 
between scalar and pseudoscalar fields. With the help of the advanced technology, the 
frequency-band selection around the optical frequency domain can be achieved by combining 
a set of high-quality optical elements such as prisms, dichroic mirrors and filters to shut out 
the non-interacting loo and uuj laser fields at the downstream of Figfj] before the detection 
of W3 with the sensitivity to a single photon. The single photon detection is not a difficult 
issue given by any conventional photomultipliers with the typical gain factor of ~ 10 6 with 
respect to the single photoelectron caused by the incidence of the single photon especially in 
the environment where coincidence signals synchronized with injections of short laser pulses 
in time are available for the rejection of the dark current noise of the photodevice. 

The brown, blue, and red lines indicate the achievable upper limits with 95% confidence 
when no photon in the frequency-band co' 3 < u 3 < tu' 3 is observed per single shot focus- 



leve 
ing 



3l| , whose parameters are summarized in TabfU We choose the laser parameters, which 
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heavily depend on the actual system, as general as possible by considering the anticipated 
pulse energies in the existing facilities reviewed above, where the wavelengths of production 
and inducing laser beams are assumed to be around 1 /zm, and the focusing parameters and 
pulse duration time are chosen so that the laser intensity [W/cm 2 ] at the surface of the final 
focusing device is lower than the damage threshold typically 10 13 W/cm 2 by more than three 
orders of magnitude considering the future pulse compression option with sub-ps duration. 
In addition, the fluence [J/cm 2 ] is also required to be lower than the typical damage thresh- 
old 10J/cm 2 for ns duration. The solid and dashed lines are for the short and long focal 
lengths, respectively. This figure provides a baseline to argue the single shot sensitivities. 
Increasing the shot statistics and shortening pulse duration time improves the sensitivity. 
These depend on the future development of the high-intensity laser technology. 

We note that the physical background process from the QED box diagram is totally 
negligible as discussed in |3| essentially due to smallness of the CMS energies in QPS. The 
generation of high harmonics from residual atoms is expected to be a background process by 
the atomic recombination process between the ejected electrons and the parent ion after the 
tunneling or barrier-suppression ionization by a strong external laser field. The appearance 
intensity values are expected to be 1.5 x 10 16 W/cm 2 and 4.0 x 10 16 W/cm 2 for iV 5+ and 6+ , 
respectively |25j]. The vacuum pressure around focal spot, therefore, should be maintained 
as low as possible. The vacuum pump commercially available can achieve ~ 10~ 10 Pa, where 
the expected number of atoms per (lOO/im) 3 volume can be below unity. We can estimate 
such a background process by requiring false combinations of linear polarization states of 
the initial and final photons in actual measurements. However, for simplicity, we assumed 
no background process in order to provide the ideal sensitivity curves at this stage. 



Filled areas are excluded domains by the other types of 



ing on the Axion-Like Particle(ALP) - photon coupling 15J as well as the upper limits 



aboratory experiments focus- 



from the searches for non-Newtonian forces by reinterpreting them based on the effective 
Yukawa interaction between test bodies 26]. The state-of-the-art methods to search for 
ALP at terrestrial laboratories by utilizing the two photon-axion coupling are represented 
by LSW(Light-shining-through-walls) 27|, the solar axion search CAST(CERN Axion Solar 
Telescope) Uand SUMICO (Tokyo Axion Helioscope) [28(, and ADMX (Axion Dark Matter 
eXpreiment) 30]. In LSW a laser pulse together with a static magnetic field produces ALP 
and the ALP penetrates an opaque wall thanks to the weakly coupling nature with matter, 
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and it then regenerates a photon via coupling to the same static magnetic field located over 
the wall. The solar Axion search is similar to LSW, but different as for the production part. 
In the Sun two incoherent photons may produce ALP and the long-lived ALP penetrates 
the Sun and the atmosphere in the earth, and they regenerate photons by coupling to a 
prepared static magnetic field on the earth. ADMX utilizes a microwave cavity immersed 
in a static magnetic field, and ALP passing through the cavity can resonantly convert into 
real microwave photons. 

The Bosonic enhancement is partly utilized in these axion searches where the axion decay 
is commonly induced under static magnetic field. However, the static magnetic field is not 
in a degenerate state with a narrow momentum range. Therefore, the enhancement of the 
decay is limited. The enhancement of the production rate is also limited because of the broad 
range of the CMS energy when choosing two photons for the production of the resonance 
state. We emphasize that the most different aspect of our approach is in the field theoretical 
treatment by which we can incorporate the nature of the resonance production and decay 
under the degenerate fields. This is in contrast to the classical treatment prescribed for 
the past axion searches. Moreover, the bulk static magnetic field has the limitation to 
increase the field strength compared with the recent leap of the laser energy [24|], where the 
cutting-edge laser technology is about to exceed Avogadro's number of photons per laser 
pulse (200kJ ~ 10 23 optical photons). 

CONCLUSION 

We have shown that the sensitivity to dark fields by searching for the four-wave mixing 
process of laser fields is expected to be able to reach the sub-eV mass range with the coupling 
strength as weak as that of gravity and even beyond it, if the cutting-edge laser technology is 
properly combined. Even before reaching extremely high fields, we have many opportunities 
to test the light cold Dark Matter candidate by the proposed method. This high-sensitivity 
is essentially realized by the Bosonic nature of laser fields. 

As a concluding remark, we emphasize some of features to search for the four-wave mixing 
process by comparing them with those in high-energy colliders as follows. As an example, 
let us remind of the Higgs production at LHC as a search for the heavy scalar field. 

First, resonance searches in collider experiments are based on measurement of the invari- 
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FIG. 4: Upper limits of the sensitivities in the (m,g/M) plane by searching for the four- wave 
mixing process with a single laser shot. The brown, blue and red lines indicate the achievable upper 
limits with 95% confidence level by assuming no background process is included, when no photon 
with frequency upshit is observed per single shot focusing, whose laser parameters are summarized 
in TabU based on the formula Eq, (|50p with C m f, = 1/2. Here T s = 2tt is used because any J-$ are 
on the same order as those in the axial symmetric case and the intention of this plot is not in the 
separation between scalar and pseudoscalar fields. The solid and dashed lines are the cases with 
short and long focal lengths. The excluded gray, black and yellow domains via the dark field - 
photon coupling are quoted from [15[. The upper limits from the searches for the non-Newtonian 



261 ] are 



force by reinterpreting them based on the effective Yukawa interaction between test bodies 
also shown with the black shaded curves. The limit of the gravitational coupling strength is added 
by assuming g ~ 0(a qe d) ~ 1CP 2 and M = Mp ~ 10 18 GeV as an order estimate based on the 
model 4|. 
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Laser parameters 


Brown lines 


Blue lines 


Red lines 


production energy 


200 J 


20 kJ 


2 MJ 


inducing energy 


200 J 


20 kJ 


2 MJ 


inducing spectrum width Su 


0.1 


0.1 


0.1 


pulse duration time r(= r c = T{) 


1 ns 


1 ns 


1 ns 


diameter at the final focusing mirror d 


40 cm 


80 cm 


800 cm 


intensity at the final focusing mirror 


3.2 x 10 s W/cm 2 


8.0 x 10 9 W/cm 2 


8.0 x 10 9 W/cm 2 


fluence at the final focusing mirror 


0.32 J/cm 2 


8.0 J/cm 2 


8.0 J/cm 2 


Short focal length / (mass range) 


1.5 m (< 0.22 eV) 


3 m (< 0.22 eV) 


30 m (< 0.22 eV) 


gm/M at 95% C.L. for short focal length 


2.52 x 10~ 27 


2.52 x 10~ 30 


2.52 x 10~ 33 


Long focal length / (mass range) 


10km (< 32 (j,eV) 


10km (< 65 (jjeV) 


10km (< 0.65 meV) 


gm/M at 95% C.L. for long focal length 


3.88 x 10~ 31 


7.61 x 10~ 34 


7.56 x 10~ 36 



TABLE I: Laser parameters to estimate the sensitivity to mass-coupling domains by searching for 
the four-wave mixing process in the vacuum. 



ant mass distribution of a produced resonance state. In searching for a resonance state in 
QPS, in contrast, we have the unavoidable CMS energy uncertainty originating from the 
uncertainty principle of optical waves compared with the extremely narrow resonance width 
due to the weakly coupling. We are unable to reconstruct the invariant mass distribution 
directly, though the interaction probability is still affected by the integrated effect over the 
possible CMS energy uncertainty. By observing the appearance or disappearance of the 
four-wave mixing signal, however, one can determine the order of the mass scale from the 
incident wavelength and the collision geometry, a significant difference from the conventional 
collider's approach. 

Second, in the case of Higgs at LHC, the dominant production channel is the gluon-gluon 
fusion process and the produced Higgs resonance state decays into two photons, where both 
the initial gluons and the final photons are not in degenerate states. Therefore, all fields 
should be treated incoherently, and the decay process occurs only via spontaneous processes 
in the vacuum, i.e., two photons in the final state are created from the pure vacuum state 
|0). On the other hand, in the case of photon-photon interactions under laser fields, all 
photons are annihilated into and created from the degenerate states. This situation results 
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in the interaction rate with the cubic dependence on the average number of photons included 
in the laser fields, which is in contrast to the square dependence of the number of charged 
particles in luminosity of the Fermionic particle colliders. The prime mission of the energy 
frontier of high energy physics is, of course, to produce new heavy particles, therefore, the 
realization of the high CMS energy is the most important task, while the sensitivity to 
weakly coupling fields is sacrificed. The dimensionless intensity included in luminosity is 
proportional to the square of the number of charged particles per bunch which is typically 
~ 10 11 particles due to the physical limitation by the space-charge effect. Even if we could 
collide them at most 1GHz over three years data taking period, the integrated dimensionless 
intensity reaches (10 11 ) 2 • 10 9 • 10 8 = 10 39 . This indicates that it is practically impossible 
to reach the sensitivity to cross sections with coupling including Mp 1 . In contrast, by the 
proposed method, we can expect the dimensionless intensity of (10 23 ) 3 = 10 69 even with a 
single laser shot including the Avogadro's number of photons. This manifestly shows how 
the proposed approach can be sensitive to the weakly coupling interactions. 

Therefore, in addition to the present most powerful experimental approach such as heavy 
boson searches at the high-energy Fermionic collider, the proposed coherent Bosonic collider 
with the inducing mechanism, simply speaking, four-wave mixing by focused high-intensity 
laser fields opens up a novel opportunity to bridge particle physics and cosmology in the 
so far unprobed low-mass and weakly coupling domains under the controllable laboratory 
environments. 

Appendix: Polarization Dependence of Scattering Amplitudes and Axial Asym- 
metric Factors Ts 

Given the scattering configuration illustrated in Figl3l the Lorentz invariant s-channel 
scattering amplitudes defined in Eq.(T2]) have the following basic form 



where S = abed with a,b,c } d = 1 or 2, respectively, denotes a sequence of four-photon 
polarization states and m is the mass of scalar or pseudoscalar field. 

The vertex factors in the numerator for the case of the scalar field exchange (SC) are 




(53) 
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defined as 

V [ f C =<P2,e ( i ) \\F, u Fn P uet ) >= \ < p 2 ,e?\F flv \0 >< 0|F""bi,e?° > 

V [ f° =< P^ef^F^^ >= \ < p„e{ d) \F^\0 >< 0\F^\p 3 ,e^ >, (54) 

while these for the case of the pseudoscalar exchange (PS) are given by 

VS PS =< p 2 ,e?\\F, u e^F pa \p l ,e ( f ) >= \ < P 2,e?\F^\0 >< 0\e^° Fpa \ Pu e ^ > 
V'f 5 =< p^ef\\F^°F prj \p^ef >=\< p^ef\F^ >< 0\e^F p(7 \p 3 ,e^ > .(55) 

Let us define the polarization vectors and momentum vectors for four photons in Figj3] 
as follows: 

e 4 {1) = (0,1,0), (56) 

(2) (2) 

e\' = (— cos-i?, 0, sini9), e 2 = (— cos^, 0, — sin??), 
e 3 = (— cos # 3 , 0, sin^), e\ = (— cos6*4, 0, — sin 6*4), 
Pi = (co sm{}, 0, uj cos??; u>), p 2 = (— usmfi, 0,a;cos$; u), 
p 3 = (co 3 sin 6*3, 0,073 cos 6*3; cj 3 ), p 4 = (— U4 sin #4, 0, U4 cos 64, U4). 



(57) 



Based on these vectors, let us summarize basic relations between momenta and polarization 
vectors with photon labels i = 1, 2, 3, 4 as follows 

(pief) = (58) 

for the coplanar condition where the plane determined by p\ and p 2 is the same as that of 
p 3 and p 4 , 

(e«e«)=l, and (e«ef ) = 0, (59) 

for any pair i, j, and 

(ef'ef') =1 for i = j, 
(efef) = cos2tf, (efef) = cos(# 3 + B A ) = cos6 + , 
(efef) = cos(tf - 6 3 ), (e? ] ef) = cos(tf - 4 ), 
(eS 2) ei 2) ) = cos(tf + 4 ), (4 2) ef) = cos(tf + 9 3 ). (60) 
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We then introduce a clock-wise rotation of the P3-P4 plane from the pi-p% plane defined 
on the x — z plane by the azimuthal angle ip varying from to 2tt around the z-axis in order 
to discuss the axial symmetry of the scattering process, when polarization vectors are fixed 
in an experiment. The rotated vectors are defined as 



Pz{ ( p) — ( w 3 sin $3 cos V 9 ) ~ u 3 sin $3 sin u 3 cos ^3! ^3) 
p^if) = (— (^sin^cos^, ^sin^sin^a^cos^u^). 



(61) 



and these result in 

(p 3 (v?)p4(v?)) = u^c^cos # + — 1) = w 2 (cos2$ — 1) (62) 

where the last equation is obtained from (pi + p 2 ) 2 = (p-3 +P4) 2 . 

With vectors defined above, the vertex factors for the scalar case are expressed as 

VS S ° = (PiP 2 )(ei Q) e^)-( Pie r ) )(p 2 ef ) ) 

V£ 1SC = {pMp^)){efe^ ) - (p 3 ^)e{ c) )(p^)e { 3 \ (63) 
and these for the pseudoscalar case are expressed as 

y [l]PS _ _ f fiupa (a) (b) 

V% ]PS = ~e^p,{ V )Mv) P eM}- (64) 

We are now ready to estimate the factor included in the partially integrated cross 
section in Eq. ffTTj) . First, we estimate S = 1122 for the scalar exchange. From the first of 
Eq.( )63|) . we obtain 

,[1]SC 



With 



we get 



V[ L f = u 2 (cos 2$ - 1) = Ku\ (65) 

(p3(^p)ef^ = — (sin 63 cos 64 cos (p + cos 6* 3 cos 64) 

(pA{*p)ei^ = 0J4 (sin 6*4 cos 6*3 cos tp + cos 84 cos 6*3), (66) 

(p3{<p)ef^ [p4(ip)ef^ = — W3W 4 {(1 + cos 2 </>) sin #3 sin # 4 cos #3 cos #4 

+ (sin 2 9 3 cos 2 64 + sin 2 64 cos 2 6*3) cos <£>}, (67) 
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and then the second vertex factor in the second of Eq. (|63|) is expressed as 

V22 = ^^{cos # + (cos 9 + — 1) + (1 + cos <£>) sin 9 3 sin # 4 cos 9 3 cos #4 

+ (sin 2 # 3 cos 2 64 + sin 2 #4 cos 2 #3) cos <^}. 

For (p — 0, this coincides with | V}i | via the relation 

V 2 2 (v 9 = 0) = ^^{cos 6* + (cos # + — 1) + sin # + } 

= W3W 4 (1 — cos# + ) = w 2 (l — cos 2$) = —/Co; 2 . 

For a small $ we take the following approximation: sin6* 3 ~ \[ 7 RjQ and sin 64 
from Eq.(l37"l) and ( 1391) . and this results in cos# + as 

cos# + ~ 1 -0 2 (1 + K), 



(68) 



(69) 



(70) 



with 



We then approximate Eq.f l68|) as 



V [2]SC 
v 22 



W 3 U;4l? 2 {cOS 2 </3 + (2cOS(y9 — 1)7£} = CJ3(x>4'{? 2 .F(y9). 



This is consistent with the approximation of Eq. (169]) for a small $ 

VM SC (<fi = 0) ~ WgW^l + £)• 



(71) 



(72) 



(73) 



By taking the square of the factorized second vertex factor, we then naturally define the 
factor JFs for the scalar case 



sc 
1122 







^ F 2 (ip)dip ~ 2vr ( | + 3TZ 2 - K ) . 



(74) 



Second, let us estimate S=1212 for the pseudoscalar field exchange as follows. From the 
first of Eq. ( 1641) with the vector definitions above, we obtain the first vertex factor as 



V Y2 



-T)-, Tin pMPVpW 



cos?9) + e pypz {- sin ■&) 



P2p 



(p w e 0ypx + p lz e zypx ) costf + (p w e 0ypz + p lx e xypz ) sin0 



P2 P [ (-coe 0ypx + co cos ®e zypx \ cos d + (-tue 0ypz + uj sin § e xypz ) sin 1? 



-we ^^ + ^ cos 



P20 cos 



0+ (- 



-we 



p 2a; + ^ sin ■&e xy0z p 2z ) sin 



^— cos + cos 0) cos + (— sin •§ — sin 0) sin *& 



-2u; 2 sin 2 0. 



(75) 
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We also get the second vertex factor from the second of Eq. (|64|) with the vector definitions 
above as follows 

= - (-e™ px cos 6 A - e pypz sin 9 A ) p 3fl p Ap 

= ' [e Qypx cos 9 A + e 0ypz sin 4 ) p 30 + (e zypx cos 6 A p 3z + e xypz sin 6 A p 3x )] p Ap 
= -cu 3 (e yzx p Az cos 6 A + e yxz p Ax sin 4 ) + 4 (e zy0x cos 0^ + e^Pa.) Pao 
= — lo 3 uj a (cos 2 6 A + sin 2 6 A cos ipj + (— cos # 4 cos 6* 3 + sin 6* 4 sin 6 3 cos y) 
= —u 3 u A cos9 A (— cos 6*3 + cos 6*4) + sin #4 (sin $4 + sin# 3 ) cosy . (76) 
For (p = 0, we find 

vjf ^(v 9 = 0) = — u 3 u A (— cos 6*3 cos 6*4 + cos 2 6*4 + sin 2 #4 + sin 9 3 sin #4) 

= -u 3 u A (1 - cos 0+) = -cu 2 (1 - cos 20) = /Cw 2 . (77) 

If we use the same approximations as the scalar case, the second vertex factor is approx- 
imated as 

V{ 2]PS ~ -w 2 u; 2 tf 2 {ft + {TT l + 1) cosy} = -u%ul& 2 G(<p), (78) 

with 

1z = -{n-Tr x ). (79) 

This is consistent with the approximation of Eq. (I77|) for a small 

y{ 2]P5 (y = 0) ~ -w 3 w 4 tf 2 (l + K). (80) 

Again by taking the square of the second vertex factor, we then naturally define the factor 
Fs f° r the pseudoscalar case 

7f 2 f 2 = P&{<p)<kp ~ 2n{n 2 + \{K- 1 + l) 2 }. (81) 
Jo 2 

Let us confirm relations for the case of (p = as follows: The ratio of the invariant 
amplitude of the pseudoscalar case to the scalar 

Mfj 2 ( V = 0) _ Vr s Vr 5 _ -2sin 2 ^ 2 ■ /G, 2 

MfUv = 0) V$ SC V® SC /Co; 2 --/Co; 2 18 j 

for a low-mass case with a small 0. The other non-vanishing invariant amplitudes are 
limited to S = 1111, 2222, 1122, 2211 for the scalar exchange and S = 1212, 1221, 2121, 2112 
for the pseudoscalar case. These relations can be confirmed by repeating routine calculations 
performed above. 
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